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Abstract 

An exact model is proposed for a gray, isotropically scattering planetary atmosphere 
in radiative equilibrium. The slab is illuminated on one side by a collimated beam 
and is bounded on the other side by an emitting and partially reflecting ground. We 
provide expressions for the incident and reflected fluxes on both boundary surfaces, 
as well as the temperature of the ground and the temperature distribution in the 
atmosphere, assuming the latter to be in local thermodynamic equilibrium. Tables 
and curves of the temperature distribution are included for various values of the 
optical thickness. Finally, semi- infinite atmospheres illuminated from the outside or 
by sources at infinity will be dealt with. 

Key words: Planetary atmospheres; Radiative transfer; Isotropic scattering; Gray 
opacity; Radiative equilibrium; Temperature distribution. 



1 Introduction 



The modern phase in modeling radiation transfer through planetary atmo- 
spheres started with the seminal work of Humphreys [1] , and the contempo- 
rary elaborate attempt of Gold [2] , followed by the critical reappraisal of both 
attempts by Emden [3] . Of these works, Emden's is directly based on previous 
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work by Schuster [4] and Schwarzschild [5], which, with Schwarzschild's second 
paper [6], laid the foundation for the theory of radiation transfer. While those 
three papers may doubtlessly be regarded as ground-breaking, they necessar- 
ily remained rudimentary from our present vantage point in radiative transfer 
theory. It was not until the next decade that Milne [7] would deal with ra- 
diative energy transfer through planetary atmospheres in a mathematically 
rigorous manner. His important developments would not be superseded be- 
fore the second half of the twentieth century, with the contributions by King 
[8], Ostriker [9], Wildt [10], Schultis and Kaper [11], Stibbs [12] and Barkstrom 
[13], to name but those authors who attempted to model exactly the transfer 
through gray atmospheres. 

By gray atmosphere we understand, as is customary, an atmosphere whose 
opacity does not depend on the frequency of the electromagnetic radiation in- 
volved. A semi-gray atmosphere is, on the other hand, one in which the opacity 
assumes just two different values in adjacent frequency intervals, conveniently 
distinguished as the "visible" and the "infrared" regions of the spectrum. It 
was Emden [3] who established this distinction, which, besides being more re- 
alistic, is susceptible to analytical treatment. The present paper is concerned 
with gray atmospheres in radiative equilibrium, while two further papers will 
be devoted to the case of semi-gray atmospheres. 

In defining the model features, we shall assume the usual ingredients of the 
classical theory of planetary atmospheres to hold true: a plane-parallel at- 
mosphere, subject to both illumination by a parallel beam of light impinging 
on one of its boundaries and to radiative heating through the other bound- 
ary, which is assumed contiguous to a partially reflecting and gray-emitting 
ground. Our purpose is to take the analytical developments as far as possible, 
under the assumption of a gray or semi-gray atmosphere in radiative equilib- 
rium. We thereby wish to push the line of investigations recalled above farther, 
drawing from the theory of radiation transfer as it evolved in the last three 
decades. Our goal is to develop elementary but rigorous models which may be 
used not only to validate more realistic models based on numerical resources, 
but also to foster the understanding of radiative transfer in more comprehen- 
sive models of, say, the climates of planets, where it plays an important part 
in long-term evolutions. 

The paper may be outlined as follows. Section 2 introduces the standard prob- 
lem of planetary atmospheres. In Section 3, the fundamental integral equation 
underlying this problem is solved. Then follows the calculation of radiative flux 
densities at both boundary planes (Section 4), which enables the evaluation 
of the ground temperature (Section 5) and the atmospheric temperature (Sec- 
tion 6) of the planet. Next, in Section 7, we calculate the radiance of radiation 
leaving the atmosphere, as it represents the fundamental quantity in solving 
the inverse problem of the theory of planetary atmospheres. In the subsequent 
section, we deal with semi-infinite atmospheres heated either from the outside 
or from within by sources at infinity. Finally, we conclude the work with a 
summary, addressing at the same time some issues to be taken into account 
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in forthcoming work. 



2 The problem stated 

We consider a plane-parallel, isotropically scattering atmosphere subject to 
the following two conditions: 

(i) its opacity is gray, i.e., independent of frequency, and 

(ii) radiative transfer throughout the atmosphere is to take place in local 
radiative equilibrium, which means that the flux density F integrated over 
frequency is independent of position. 

To these defining physical assumptions we add the two following boundary 
conditions: 

(Hi) the atmosphere's upper surface is illuminated by a collimated solar beam 
(corresponding to a sun at infinity), carrying a steady frequency- integrated 
flux Fq through a unit surface normal to the beam, while 
(iv) the bottom boundary of the atmosphere is in contact with a gray, isotrop- 
ically emitting and reflecting ground. 

Note, however, that in Section 8, devoted to semi-infinite atmospheres, the 
boundary condition (iv) will be replaced by an appropriate condition at infin- 
ity (v). 

Turning to the basic equation of the model, let us begin by first introducing 
the frequency-integrated radiance field appropriate to our model assumptions. 
This function I(r,fi,(p) depends on a positional and two angular arguments. 
The optical depth variable r covers the range [0, r 6 ], where r h is the optical 
thickness of the atmosphere. The top of the atmosphere is at r = and the 
bottom (i.e., the ground) at r = tj,. t and are frequency-independent, /i 
stands for the direction cosine of the angle of incidence (zenith angle), referred 
to the outward normal at r = 0, and ip is the azimuth angle with respect to 
a chosen tangent direction in the plane r = 0. The radiance field obeys the 
following radiative transfer equation 



wherein the source function S is composed of the thermal emission and the 
radiation scattered isotropically by the elementary volumes of the atmospheric 
gas. The thermal emission may be written in the form [1 — a(r)]23*(r), where 
a(r) denotes the albedo for single scattering, which is independent of frequency 
by assumption, while 03* (r) is an unknown function of r, deriving from the 
Planck function if the atmosphere is in local thermodynamic equilibrium. The 
scattering source function is expressed as a(r)J(r), where J(r) is the mean 
radiance, defined by Eq. (6) below. Hence 






(2) 
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The boundary condition on the (frequency-integrated) radiance at the upper 
surface r = is: 



7(0, /x, y?) = F 5(fi + fi )5(<p - <p ) (-1 < fi < < fi < 1), (3) 

where the angular variables (/xq, <£>o) define the direction of incident solar ra- 
diation, while S denotes the Dirac distribution at 0. /x is positive because 
it refers to the zenith angle of the source (a point-like sun). At the bottom 
surface r = r& we require 

I(r b ,fi, V ) = -F\r b ) (0</i<l), (4) 

71 

the upward flux density at the bottom level r = r& being given by 

F\r b )=F}+r s F^r b ). (5) 

The upward and downward fluxes are defined by Eqs. (9) and (10) below, 
respectively. Fj refers to the flux radiated by the ground, to be calculated in 
Sec. 5, and r s is the Lambert flux reflectance of the ground (with its usual 
meaning, see e.g. [14]), which we assume to be given. 
We next recall the definitions of mean radiance and net flux density 

I r2ir r+l 

J ^ = 4nh l-i / ( r ' /i ' V9 ) d/idv? ' ( 6 ) 

r2n r+l 

i? ( r ) = y J I ( T >i i > < p)i j di id( p> ( 7 ) 

all quantities being integrated over frequency. Equation (7) may be rewritten 
in the form 

F = F^t) - F 1 (t), (8) 
in which the upward and downward flux densities are respectively given by 

^ T (r) = r /' /(r, /i, y?)M"dy? > 0, (9) 
Jo Jo 

Fi (.r) = -J Q j^I^ti^fjdfidfp^O. (10) 

We recall that the condition of radiative equilibrium of a gray atmosphere 
means that the net flux density F(r) is a constant F, which appears on the 
left-hand side of Eq. (8). 

The formal solution of the radiative transfer equation (1) under the boundary 
conditions (3) and (4) reads 

1 f T 

I(t,h,<p) =F 5(i2 + fi )S((p - (fo) exp(r/» / S(t') exp[(r - r^/^dr' 

u Jo 
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for — 1 < /i < 0, and 
I(T,n,<p) = V(r 6 )exp[-(r 6 -r)//x] + - P S{r') exp[-(r' - r)/ M ]dr' (12) 

7T /i Jt 

for < /x < 1. We may thus derive the following expressions for the mean 
radiance and the two flux densities previously defined, in terms of the source 
function: 

J(r) = J,{r) + \ £ El (\r- r' \) S (r')dr' , (13) 
F\r) = Fq (r ) + 2n£ b E 2 (t' - T)S(r')dr', (14) 

F 1 (t) = Fi(r) + 2vr f E 2 (r - r / )-5(r / )dr / , (15) 
Jo 

where Jo( r )> ^o( T ) an d -^o( r ) denote the corresponding contributions from 
the two boundary planes: 

J (r) = ^ exp(-T//! ) + ^^(r,)^^ - r), (16) 

F T (T) = 2FT(r b )£ 3 (r 6 -T) , Fj(r) = exp(-r/// ). (17) 

The functions E n appearing on the right-hand sides of Eqs. (13)-(17) stand 
for the exponential integral functions of order n > 1, defined by 

E n (r)= / 1 exp(-r//i)/i n - 2 d/i. (18) 
Jo 

The main purpose of the present article is to calculate the incident and re- 
flected fluxes at both boundary planes, as well as the temperature of the 
ground and the temperature profile of the atmosphere, in terms of the inci- 
dent flux F and the constant flux density F. Note that F is truly a constant 
in our model, whereas F depends on the position /z of the sun, as do all other 
quantities describing the radiation field. Beginning in the next section, this 
dependence will be made explicit. 



3 The fundamental integral equation solved 

The condition of radiative equilibrium entails that the atmosphere absorbs 
locally as much radiative energy as it emits, and thus 

J(t,//o) = S(t,// ) = «*(t,//o) (19) 

for a gray atmosphere. The first equality is equivalent to the condition F(t, /i ) 
= F(fi ), as may be seen from the integration of both members of the radiative 
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transfer equation (1) over n between -1 and +1. Replacing J by S on the left- 
hand side of Eq. (13), we obtain the following integral equation for the source 
function: 

S(t,im>) = J (t,^) + 1 - H E^t-t'^t'^vW. (20) 

I Jo 

This equation shows that the frequency-integrated source function solves an 
integral equation of the Schwarzschild-Milne type with albedo equal to 1. 
Physically, this is due to the fact that a gray atmosphere in radiative equi- 
librium is conservative for the frequency-integrated quantities describing the 
radiation field: there is neither loss nor gain of energy by absorption and 
emission, which is what the first equality in Eq. (19) is telling us. This fact 
by no means implies that the atmosphere's albedo for single scattering is 
everywhere equal to 1. Quite the contrary is the case when the second equal- 
ity in (19) is imposed, which may be seen by inserting the first equality of 
Eq. (19) into the definition (11) for the integrated source function. One gets 
[1 — ci(t)]S(t) = [1 — a(r)]5*(r), leading to the second equality in (19) for 
layers r in which a(r) ^ 1. Once the term 1 — a(r) drops out, the albedo for 
single scattering of the atmosphere plays no role in the model, which means 
that it can take on any value, and in particular be spatially variable. 
In view of Eq. (16) and by the linearity of equation (20), its solution is found 
to be 

S(r, !*>) = -[ ^B(r, &>) + F\r b , fi Mn - r) ] , (21) 

71 4 

where the functions B and £ solve the following integral equations: 

£(T,//) = exp(-T///) + \ p E^r - t'\)B(t> \pt)dr\ (22) 

2 Jo 

Mr) = \e 2 {t) + \ £ Ei(\t — r'\)Ur')dr'. (23) 

The functions B and £ are classical auxiliary functions of transfer theory, 
that have to be calculated for a slab of single scattering albedo equal to unity 
and an arbitrary optical thickness r^. In what follows, their dependence on 
these two parameters will be implied. Thus, we shall write B(t, /i) instead of 
B(l, r b , r, /i), £o( r ) instead of £o(l, r b, r ), an d likewise for all the other auxil- 
iary functions we shall introduce. As we just saw, the fact that these auxiliary 
functions are to be calculated for a conservative slab (i.e., of unit albedo) is 
without relevance to the real value of the atmosphere's albedo, which may at 
any point be arbitrary. 

The function B was first introduced by Ambartsumian [15] in his attempt 
to solve the problem of diffuse reflection from a semi-infinite atmosphere. It 
was studied in detail by Busbridge [16] and coincides with van de Hulst's 
point-direction gain, which clarifies its physical meaning [17]. An analytical 
expression for it in the case of a conservative atmosphere is given in the Ap- 
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pendix. The function B satisfies the following relation referring to a finite slab 
of thickness r b 

B{t,h) = B(r h - r, -//)exp(-r b ///), (24) 
and it reduces to the well-known X- and F-functions on the two boundary 
layers: 

B$,n)=X(n) , B{r b ,n) = Y( i i). (25) 
We refer to Chandrasekhar [18] or Busbridge [16] for a survey of the principal 
properties of the functions X and Y, along with their moments of order n > 



a n = f X(fi)fi n dfi , /3 n = f Y(fi)fi n dfi. (26) 
Jo Jo 

As to the function £o ; it may be seen from Eqs. (22)-(23) that it is proportional 
to the angular moment of order of the function B, i.e., 

Ur) = \ f B{r^W. (27) 
2 Jo 

Equation (23) also shows that it can be interpreted as the escape probability 
function through the top surface of the atmosphere [17]. It follows that 

Ur)+Un-r) = l, (28) 

which means that a photon, emitted at depth r within a conservative finite 
slab, will eventually escape through either boundary surface. Setting r = 
in this relation and noting that the boundary values of the function £o are, 
according to Eqs. (25)-(27), given by 

£o(0) = ^ao , Un) = \(3o, (29) 

one may readily deduce that the moments of order of the functions X and 
Y fulfill the relation a + f3 = 2 in conservative media. 

We refer to [17] for graphs and tables of the function g = 2£ , and to [19] for 
the analytical calculation of the function £ 0) including very accurate tables of 
this function. In the Appendix, we reproduce its expression for the conservative 
case. 

In the remainder, we will require the finite Laplace transform, at (fj, > 0), 
of the function r — > B(r,fi ). Its expression 



Jo 



Tb 

J B(r,/i )exp(-r/ / u)dr = 4/i/i -R(/i, /i ) (30) 

o 



can be found in [16], p. 90. The reflection function R(/i, fj, ) appearing on the 
right-hand side is the usual function describing the radiation reflected by a 
slab illuminated on its upper surface: see, e.g., Sec. 4.1.2 of van de Hulst [17]. 
The latter reference produces the following expression for a conservative slab 
(p. 194): 

%/l) = 4 ^ • (31) 
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Integrating both members of (30) between and 1 with respect to /io, H 
and both variables in turn leads to the following exact relations involving the 
functions B and £ ; 

P &(r) exp(-T/»dr = /li^/i), (32) 
J o 

\ l Tb B(r,fi )E 2 (r)dr = fi RM, (33) 
I Jo 

rto(T)E 2 (T)dT = ±R 11 , (34) 
Jo 2 

where i?i(/i) and i? n denote the moment and bimoment of order 1 of the 

reflection function, respectively, viz. 

RM = 2 f 1 R{n, //)//d// , R u = 2 C RMfxdfx. (35) 
Jo Jo 

The calculation of these coefficients from expression (31) of the reflection func- 
tion is well known: see, e.g., Display 9.1 on pp. 194-195 of [17]. The result is, 
in the conservative case, 

Ri(fi) = l-^Po(X + Y)(n) , R n = l-/3o(ai + /3i). (36) 

Now, by virtue of the properties (24) and (28) of the functions B and £o, the 
following relations may be derived from Eqs. (32)-(34): 

rn 

/ ^o(^-^)exp(-r//i)dr = /i[T 1 (/i)-exp(-r 6 //i)], (37) 
Jo 

rn 

/ B(r b - T,Ho)E 2 {T)d.T = 7i(A*o) - exp(-r 6 //x )], (38) 
Jo 



l r* 

2 



r £o(r 6 - r)T 2 (r)dr = - [T n - 2£ 3 (r 6 )]. (39) 
Jo 2 

Here Ti(/x) and Tn are the moment and bimoment of order 1 of the transmis- 
sion function of the atmosphere, whose expression is 

[see Sec. 4.1.2 and Display 9.1, p. 194-195, of [17]]. We thus have 

T.ifi) = 2 /' 1 T( / x, / /) / /d// , T n = 2 f' TMfidfi, (41) 

JO JO 

and from Display 9.1 of [17] 

T 1 ( l i) = ±p (X + Y)(ri , T 11 = A)(ai + A)- (42) 
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We see that indeed Ri(fi) + Ti(//) = 1 and Rn + Tu = 1, as expected for a 
conservative atmosphere. 

The solution (21) is, at this stage, merely formal, as it depends on the yet 
unknown flux density F^(r b ,fi ). In the following section, we shall calculate 
this flux by replacing the source function, given by expression (21), into Eqs. 
(14) and (15) for r = and r = r b) respectively. 



4 Calculation of the flux density through the boundary planes 

According to Eqs. (14)- (15) and (17), the upward flux at r = and the 
downward flux at r = r b are respectively given by 

F\0,fi ) = 2F\r b ,fi )E 3 (r b ) + 2n P E 2 (t)S(t, /x )dr, (43) 

Jo 

F l (r b , fi ) = fi F exp(-r b / fi ) + 2tt / " E 2 (r b - t)S(t, // )dr. (44) 

Jo 

Substituting expression (21) for S(r,fi ) in these integrals and taking into 
account Eqs. (33)-(34) and (38)-(39), one obtains: 

F T (0, fio) = RiMfMiFo + T u F\n, fi ), (45) 
F i (r b ,fi )=T l (fi )fi F + R ll F\r b ^ ). (46) 

The flux densities at the top of the atmosphere have values that follow im- 
mediately from the boundary condition (3) and flux constraint (8) at r = 0: 

i^(0,/i ) = fi F , F T (0,/i ) = fi F + F(fi ). (47) 
We infer the two fluxes at the bottom boundary plane from Eqs. (45)-(47) 

F\r b ,fi ) = -L^FoTM + F(/i )], (48) 
1 n 

F l (r b ,fx ) = J-^FoTl^o) + RuFfa)]. (49) 
1 n 

Note that T n ^ if the atmosphere is finite {r b < oo), which is what we are 
supposing in this section. 

These relations express the four boundary fluxes in terms of the incident flux 
Fq and the flux constant F(fi ). They have the virtue of simplicity when it 
comes to calculating the surface fluxes when at its top the atmosphere receives 
in a unit of time as much energy as it radiates [F(/i ) = 0]. At the same time, 
their physical interpretation is not as easily grasped as when the four boundary 
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fluxes are expressed in terms of the incident flux Fq and the flux emitted by 
the ground, F](/j, ). To see this, we express the constant flux F(fi ) in terms 
of the fluxes F and Fj(/x )- But first we note that relation (8), for r = r b , 
and relation (5) together imply 

F^o) = F}(fi ) - (1 - r s )FVfc,/io), (50) 

in which we now may substitute the expression (49), thus giving, after some 
rearrangements 

F(fio) = T*FJ(hq) - (1 - r s )T(n )n F . (51) 

T(hq) and T* denote respectively the flux transmittance and the spherical 
transmittance of the atmosphere, whose definitions will be recalled further 
on. 

Replacing this expression for F(/i ) in Eqs. (47)- (49), we find the following 
relations: 



FT(0, fi ) = AM^F + T*F](no), (52) 
F\r b , f i )=r s T( f i ) f i F + (T*/T 11 )F}( f i ), (53) 
F\r b) ^) =T(fM>)tM)Fo + T*(R 11 /T 11 )F}( f i ), (54) 

in which appear the flux reflectance (or plane albedo) A(fio) and the flux 
transmittance (or transmission) T(fi ) of the atmosphere, whose expressions 
for a conservative atmosphere are [20,17] 

A(^ ) = 1 - (1 - r.) ri(/io) = RM + r Tl \ Ti(/x ), (55) 

™ = r^- (56) 

The associated spherical quantities 

A* = 2 f 1 A(/i )/i d/io , T* — 2 f 1 T(/i )/i d/io (57) 
Jo Jo 

are thus given by 



- 4 '- 1 -( 1 - r -»T^k^" +r T^k' < 58 » 
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Note also that the reflectances and transmittances are not independent of each 
other, as they verify 



A(no) + (1 - r a )T{fio) = 1 , A* + (1 - r s )T* = 1 (60) 

in the case of a conservative atmosphere. Equation (51) clearly relates inso- 
lation of the planet, emission of its surface, and constant flux density F(/i ) 
throughout its atmosphere. Notice that the latter is proportional to the co- 
efficient f3o, which approaches zero when the optical thickness r b tends to in- 
finity. The constant flux density -F(/xo) necessarily vanishes for a semi-infinite 
atmosphere irradiated from the outside, which is to be formally confirmed 
in Subsection 8.2. Equations (52)-(54) extend the usual relations for atmo- 
spheres overlying an emissive ground. It will be noted that (53) reduces to 
F^(r b ,fi ) = Fj(fi ) when the ground is totally absorbing (r s = 0). 
We introduce now what will be called the atmosphere's "grayness factor", 
defining it as the quotient between the flux density arriving at the ground and 
that reaching the top of the atmosphere 

G{no) = — -— — . (61) 

It may be rewritten in two different ways, depending on which expression is 
laid down, (49) or (54). The first one, 

G W = ^ + (62) 

-ill -ill ^0-^0 

enables us to calculate the factor readily when the atmosphere receives as 
much energy as it radiates in a unit of time [F(/i ) = 0], in which case 

in 2 «i+Pi 

This factor increases with /i , as of course it should. It does not vanish for 
grazing incidence of the solar rays (/xq = 0), but remains always smaller than 
0.5. For normal incidence (/i = 1), whatever t&, it is greater than unity, 
being enclosed within the limits 1 (r b = 0) and approximately 1.259, the value 
corresponding to r b = oo (Subsection 8.2). For intermediate values of /i , G(/i ) 
may be smaller or greater than unity, depending on the values of r b , as may be 
gleaned from published tables for the functions X, Y and their first moments. 
Such tables may be found, for example, in [21]. The grayness factor is greatest 
for semi-infinite atmospheres irradiated from the zenith (r b = oo, fi = 1). 
The other expression for the grayness factor is obtained if in (61) we replace 
the flux reaching the ground by its expression (54), giving 

GM=TM+T*|li^^. (64) 

ill Ho-Fo 
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We see that this coefficient generalizes the transmission coefficient of the atmo- 
sphere, with which it is identical when the planet's surface does not radiate 
thermally. Moreover, the grayness factor provides an interesting alternative 
to the relation (51) linking the thermal emission of the ground to the flux 
constant F(/x ), as may be seen at once from Eqs. (50) and (61): 

F(/i ) = Fj(/i ) - (1 - r s )G(fi )fioF . (65) 

So as to be able to describe the global characteristics of the atmosphere, i.e., 
the ones at the planetary scale with a half-sphere irradiated by the sun, it 
would seem judicious to introduce a spherical grayness factor in the same way 
as the definition of the spherical transmission. The flux density of the radiation 
reaching the ground of the illuminated hemisphere is 2ttR 2 Jq 1 F^u, /! )d/io, 
where R is the radius of the sphere. Dividing this flux by the solar flux F irR 2 
intercepted by the sphere, we obtain the following definition of the spherical 
grayness factor: 

G* = ^- I' Fl(n, /i )d/i = 2 C G^oVocW (66) 

-To Jo Jo 

This average coefficient fulfills a number of interesting relations, obtained by 
integrating Eqs. (62) and (64) over /i between and I. Thus arises the globally 
averaged flux constant 

F = \[ F ^)^o (67) 
and likewise for the mean thermal flux emitted by the ground 

F} = \ £ F*<jm>)imAiM>. (68) 

The factor 1/2 in these expressions represents the ratio between the irradiated 
surface and the total surface of the atmosphere. 

The formulas for the spherical grayness factor deriving from Eqs. (62) and 
(64) read: 

(T^ + ^^Hl + Afi], (69) 

-ill *0 -ill ro 

and by integrating (65) over /i one gets the relation 

F = F]-(l-r s )G*^. (70) 

The first equation in (69) shows that the spherical grayness factor is equal to 
1 when the atmosphere is in global radiative equilibrium, that is, when the 
global average flux F vanishes. The second equality in (69) shows that the 
spherical grayness factor extends the concept of spherical transmissivity to 
the case of heating from below, both expressions being identical when there is 
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no such heating by thermal radiation from the ground. Finally, an interesting 
variant of (70) arises from the elimination of the ratio Rn/Tu between the 
two equations in (69), giving 

(G* - l)T*F} = (G* - T*)F. (71) 

We shall take advantage of this result in the next section, devoted to calcu- 
lating the temperature of the planet's surface. 



5 The ground temperature calculated 

To derive the temperature T s (/i ) of the ground from the flux Fj(/^o) it radi- 
ates, we have to adopt an emission law for the ground. The simplest one is 
the gray law with (flux) emittance e s , viz. 

e s aT^ )=F}(fi ), (72) 

where a denotes the Stefan-Boltzmann constant. The flux Fj(/i ) radiated by 
the ground may be evaluated from Eq. (65), thus enabling us to express the 
temperature of the ground in terms of the grayness factor: 

e s aT*(pi ) = (1 - r s )G(fi )fi Fo + F(pt ). (73) 

This simplifies when the flux constant F(/i ) vanishes, and the equilibrium 
relation e s + r s = 1 is assumed to hold true, becoming 

aT s Vo) = GMfjioFo. (74) 

Note that in this particular case the surface temperature is independent of the 
emitting and reflecting properties of the ground. 

The temperature thus defined depends on the sun's elevation, which is to say 
that it assumes different values according to the location on the illuminated 
hemisphere. One may also define a global-mean surface temperature T s of a 
planet with radius R in rapid rotation by equating the flux emitted by its 
surface, 47iR 2 e s aT^, to the flux irradiated by the surface of the illuminated 
hemisphere, 2irR 2 Jq F} (fx ) fx dfi = 2irR 2 e s a Jq T s 4 (/i )// d// - Therefore, 

1 r 1 

aT ° = 2Jo aT s(Vo)Vo d Vo- ( 75 ) 

It is thus easy to deduce from (73) and the first equality in (69) the following 
expression for the mean thermal flux emitted by the planet's surface: 

e s oTt = (1 - r s )^G* + F = (1 - r s )^f + (76) 



13 



For an atmosphere in global radiative equilibrium (F = 0), bounded by a 
surface for which e s + r s — 1, the surface temperature is simply given by the 
relation 

oTt = ^, (77) 

which shows that it depends only on the planet's distance from the sun. 
As an illustrative numerical instance, let us determine the main global char- 
acteristics of the terrestrial atmosphere, assumed to be gray and in radiative 
equilibrium. It is known that the solar constant F has a value of 1370 W/m 2 , 
and that the spherical albedo A* is close to 0.3. Adopting the mean value 
r s = 0.095 for the gray reflectivity of the Earth's surface (r s = 0.15 for visible 
light and r s = 0.04 for infrared radiation, see [22]), the value of the coeffi- 
cient of spherical transmission may be deduced from the stated value of the 
spherical albedo: T* = 0.77. We note in passing that the knowledge of r s and 
either A* or T* allows us to infer the optical thickness r b of the atmosphere: 
for that to be possible, it suffices to solve one of the equations A*(r&, r s ) = 0.3 
or T*(r{,,r s ) = 0.77 by bringing to bear (58) or (59). One finds r b = 0.379 
as the optical thickness of the Earth's atmosphere when it is supposed to be 
both gray and in radiative equilibrium. 

The mean surface temperature T s is observed to have a value around 288 
K [14,22], and the coefficient of emission of the surface is found to be e s = 
l — r s = 0.905. One may therefore deduce from the second equality in (76) the 
value of the mean flux constant F = 33 W/m 2 , then from the first equality 
the value of the spherical grayness factor: G* = 1.03. The ratio F/F is of the 
order of 2%, which is close to a state of global radiative equilibrium. Let us 
remind the reader that strict global radiative equilibrium is incompatible with 
the equilibrium relation e s + r s = 1 in our model, for it would imply relation 
(77), which is not satisfied by the chosen values for F and T s . 



6 The atmospheric temperature calculated 

If now we assume the atmosphere to be in local thermodynamic equilibrium, its 
temperature structure may be determined by substituting the expression (48) 
for the upward flux into the solution (21) for the source function S(r,fi ) = 
*B*(r, /i ) = (a/n)T 4 (r, fi ). Doing so, we obtain 

oT (r, no) = —B(t, no) + £ Q (T b - r). (78) 

This expression simplifies when F(n ) = 0. Dividing both members by F , 
defining the temperature T by <tTq = F , and recalling that Ti^^/Tn = 
G(/i ), we may express the temperature distribution within the atmosphere in 
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the following dimensionless form: 



T(t, Ho) 



1/4 



(79) 



The surface values of the atmospheric temperature may be calculated thanks 
to the surface values (25) and (29) of the functions B and £o- 111 this way, we 
get 



T(0,Ho) 
To 

T(n, Ho) 
To 



1 1 i l ' A 

-X(no) + -PoHoG(ho) 



-£ Y M + ^oHoG(ho) 



1/4 



(80) 



(81) 



In particular, with the help of Eq. (74), one can write the temperature at the 
bottom of the atmosphere in the form 



T 4 (n,fio) = -j-Y(no) + -a T*(fio) , 



(82) 



which clearly shows that it differs from the surface temperature T s (/i ). This 
discontinuity extends the well-known classical discontinuity occurring when 
F = o-T 4 = [14,19]. Since Y(ho) > and (l/2)a < 1, the temperature 
of the atmosphere at the ground level may exceed or fall short of the ground 
temperature T s (ho), depending on the specific values of Ho an d T b- Likewise, 
the following discussion will make clear that T(0, Ho) < T{r b , Ho) or T(0, Ho) > 
T(r fe ,/x ), depending on the value taken by the parameters /i and r b . 
The ratio T(r, /x )/T is tabulated in Table 1 for /x = 0.25. This value has 
been chosen as appropriate to a global and annual average of the instantaneous 
insolation of every point on a spherical planet revolving around the sun in a 
circular orbit. The ratio is plotted in Fig. 1 for fio = 0.1,0.25,0.5 and 1. The 
selected values of r b are 0.01, 0.1, 1 and 10 for each one of these values of 
Ho- We made use of Eq. (A4) in the Appendix to evaluate the conservative 
function £ - The present temperature profiles should of course be different 
from those of [19], inasmuch as (79) is based on the assumption F(/i ) = 0. 



The shape of the temperature curves shown in Fig. 1 depends on the angle of 
inclination of the solar rays. For /io = 0.1 and fio — 0.25, the temperature tends 
to decrease from a higher level- value in the surface layers of the atmosphere to 
a lower level-value in layers where the atmosphere is optically thick enough, 
which may be assumed to be the case when r b > 1. The temperature levels 
off roughly for r > 1 or 2, and the atmosphere becomes virtually isothermal 
at greater optical depths. The two "plateau" temperatures are given by the 
expressions (80)-(81), and depend but weakly on the optical thickness r b of 
the atmosphere. 

For greater values of //q, Ho — 0.5 and /i = 1 for instance, the temperature 
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Table 1 

Values of T(r, hq)/Tq for = 0.25 and r& = 0.01, 0.1, 1 and 10. 



n t 


T{t,ho)/T q 


n 
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T(t,ho)/T 


n ni n 
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D 7836 

U. 1 oou 
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u. / 00 ( 
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VJ. 1 000 
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VJ. 1 t OVJ 


0.005 
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0.05 
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0.006 


0.7828 




0.06 


0.7720 


0.007 


0.7826 




0.07 


0.7688 


0.008 


0.7823 




0.08 


0.7655 


0.009 


0.7821 




0.09 


0.7621 


0.01 


0.7817 




0.1 


0.7583 


1 n 

-L U 


D 7886 

U. 1 oou 


1 n 

J-U 


n 

VJ 


D 7886 

VJ. 1 oou 


n 1 

u. j. 


D 7649 

VJ. 1 \)'-±Li 




1 

_L 


D 65D1 

VJ .VJOVJ ± 


n 9 


n 7374 

VJ. 1 O I 




9 


D 6408 


n 3 

U . 


D 7147 

U. 1 11 1 






D 6399 

VJ .VJkJC/C/ 


n 4 








u.vjoyo 


0.5 


0.6824 




5 


0.6397 


0.6 


0.6717 




6 


0.6397 


0.7 


0.6636 




7 


0.6397 


0.8 


0.6575 




8 


0.6397 


0.9 


0.6530 




9 


0.6397 


1 


0.6494 




10 


0.6397 



starts off increasing as we get into the atmosphere, a fact that reflects the 
greater atmospheric heating efficiency of the sun's radiation at solar incidence 
angles closer to the surface normal. The temperature reaches a maximum 
value when /i = 0.5, but for normal incidence (/xq = 1) it does not cease to 
increase with depth, whatever the value of 75. The fact that T(t{,,//o) may 
be greater than T(0,/xo) for normal or nearly normal incidence of the solar 
rays may readily be understood by noting that these two extreme temper- 



16 





Fig. 1. T(t,u )/T versus r for fi = 0.1,0.25,0.5,1 and n = 0.01,0.1,1,10. The 
curves for r b = 1 are not visible since they nearly coincide with that for t& = 10 
over the range [0, 1]. 

atures are close to their values in a semi-infinite atmosphere as we plunge 
deep enough into the latter, say for r h > 1. Indeed, from Eqs. (108) below, 
we have T{t^ ii o )/T(0, fj, ) ~ T(oo, fi )/T(0, /i ) = (v^fyto) 1 / 4 , and this ratio 
exceeds unity whenever /i > 1/ y^- It is also under these circumstances that 
the grayness factor is necessarily greater than unity. 

To illuminate further the behavior of T(r, /x ) according to the values of /x 
and r b , we need to ascertain the temperature derivative at every point of the 
atmosphere, which calls for the previous determination of the derivatives of 
the functions r — > B(r,fio) and r — > £o( r & — T ), a s can be seen from (79). 
Fortunately, these two derivatives may be expressed analytically in terms of 
a new auxiliary function, denoted by 0, a function that has been thoroughly 
studied by Sobolev [23,20]. While the expression for the derivative of B is 
classic (see [16], p. 91), the derivative of £o results from it by a simple integra- 
tion over the angular variable /z. We shall dispense with these calculations for 
finite atmospheres and instead carry them through only in the case of semi- 
infinite atmospheres (Subsection 8.2). We point out that a thorough study of 
the temperature variation in the neighborhood of both boundary planes would 
reveal that the derivatives of the temperature are infinite at r = and r = r b , 
whereas the curves of Fig. 1 seem to suggest that they vanish, an illusion due 
to the chosen logarithmic scale for the abscissa. 

We now define the global mean temperature T(r) of a planet's atmosphere 
in fast rotation in much the same way as we defined its surface temperature, 
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that is, through the relation 



aT\r) = l -j\T\r,^W i 



o- 



(83) 



This temperature profile verifies the equation obtained by integrating both 
members of (78) over /x , viz. 



In an atmosphere in global radiative equilibrium, the average flux constant 
F vanishes and the preceding relation shows that the averaged atmospheric 
temperature is constant and equal to Fq/A. This is just the value encountered 
for the surface temperature when one assumes both the equilibrium relation 
e s + r s = 1 and F — 0. Both assumptions thus imply that the ground and 
atmosphere are, on a global average, thermally uniform. 



7 The radiance of the emerging radiation from the atmosphere 

To be able to calculate the radiance of the radiation field in the atmosphere, 
we need only substitute the expression (21) for the source function into the 
formal solution (11)-(12) of the radiative transfer equation, and then integrate 
analytically with respect to r' . Although such an integration is possible for 
finite media, it is rather tedious, due to the complexity of the functions B and 
£o in the finite case. It is much more tractable for semi-infinite atmospheres, 
for which we shall give the result in the next section. Here, we limit ourselves 
to recalling the expression for the emergent radiance, a quantity useful in the 
solution of inverse problems raised by the theory of planetary atmospheres. 
Inserting the expression (21) for the source function into relation (12) for 
r = 0, one can perform the integration with respect to r' thanks to Eqs. (30) 
and (37). We get for < n < 1: 



We recall that the reflection function R(fj,, /x ) is expressed by relation (31) in 
terms of the X- and F-functions. It is clear that we can infer the flux F^fa, fio) 
radiated from the ground by measuring the radiation leaving the atmosphere, 
provided we know Tj,, F and fiQ. Likewise, the ground temperature may be 
deduced from F\r b , /x ) with the help of Eqs. (72), (5) and (46), should the 
reflectance and emittance of the ground be known. 



(84) 



/(0,/i,^) = - \fi F R(fi,fi ) + F\r b ,fi )(l/2)(3 (X + Y)(fi) 

7T L 



(85) 
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8 Semi-infinite planetary atmospheres 



We may distinguish two categories of semi-infinite atmospheres: those that 
are heated from the outside and those with internal sources at infinity, as in 
the case of stellar atmospheres. Of course, that does not prevent atmospheres 
from having both categories of sources. In this event, the radiative quantities 
required for their description are the sum total of quantities resulting from 
each category treated separately. That we may proceed in this manner is a 
consequence of the linearity of the transfer equation under the given boundary 
conditions. 

As gray semi-infinite models are relatively simple, we shall be able to evaluate 
the radiance at each point of the atmosphere, a task we have laid aside in 
the finite case. To that end, new auxiliary functions will be required, whose 
analytical expressions we recall in what follows. 



8.1 The functions H , <p and rj 



The function H, first introduced by Hopf [24] for conservative media, has 
been well studied and it is discussed in most treatises on radiation transfer, 
beginning with Chandrasekhar's and Busbridge's monographs [18,16]. In a 
conservative atmosphere, the function H may be evaluated with the aid of the 
following exact formula, resulting from Mullikin's expression when the albedo 
approaches unity [25]: 



H(/j,) — (1 + jj) exp 



lo v(v + /i) 
It involves the functions defined for < v < 1 by 



(/i > 0). 



:s6) 



6{y) = — arctan 

7T 



7T V 



2T(v) 



*7) 



T(v) = l-±v\n(±±^). (88) 

Continuous values on [0, tt] of the arctan function are used in the definition 
(87) of the ^-function, i.e., the branch is not the principal one. The resulting 
function arctan(x/y) is usually denoted by ATAN2(x, y). With this choice, the 
function v — > 6{v) is continuous from [0, 1[ to [0, 1[, although T(v) vanishes 
once in the interval [0, 1[. 

The function 0, or resolvent function, is the other auxiliary function that has 
been widely studied, notably by Sobolev [23,20]. For conservative media, we 
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may calculate it by means of Minin's expression [26] 

0(r) = v / 3 + i/ 1 jjj^eM-rM- (r > 0), (89) 
2 Jo H [v) v 

where 

= p M+ 1(,m (0 £ " < r) ' (90) 

The function 77, on the other hand, is far less well-known than the preceding 
two functions, and only occasionally it has cropped up in the literature, those 
occasions having already been pointed out in Section 4 of [27]. We shall cal- 
culate it here for conservative media by letting a — > 1 in Eqs. (36) and (38) of 
that article, which yields for r > and /i£l: 

,(r,,)=Y( M m^tlM +s(T ^ (91) 

Here Y(/x) denotes Heaviside's function, equal to for ji < and to 1 for 
/i > 0, and 



c(/i) 



l/T(/x) if /i e] - oo,0[U]l,+oo[, 

^)T(ii) if /!£ [0,1[, (92) 
if fx = 1, 



and 

^'^U'li^-^^ (93) 

The integral on the right-hand side has to be evaluated as a Cauchy principal 
value for < /i < 1. For r = and /1 ^ 0, this integral may be carried out 
analytically, whereby we find 

™» = l ll -m ] - (94) 

Consequently, 

r,(0,n) = ±[l-Y(-n)H(-n)]. (95) 
jj, 

Equations (91)-(93) also show that the function 77 is bounded as r — > 00, with 
limit equal to \/3- Finally, these equations can be extended by continuity at 
H — 0, and we retrieve the function 

77(r,O) = 0(r) (r>0). (96) 

The functions H, <p and rj are the building stones that enable us to construct 
almost all auxiliary functions needed in solving the transfer equation in a half- 
space, the function B being one of them. Indeed, the latter is related to the 
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functions H and r\ by means of the following relationship, valid for /i > [27] 

B(r,fi)=H(fi)fir ] (r^). (97) 



8.2 Semi-infinite atmospheres irradiated from the outside 



In this model category, to which, as a particular instance, belongs the atmo- 
sphere of Venus, the third condition (Hi) of Section 2 applies unchanged: 
(Hi) At the upper surface of the atmosphere, the incident radiation is given 
by 

I(0,H,ip,fM),ipo) =F 5(n + /j, )5((p-(po) (-1 < /I < < fj, < 1), (98) 

while condition (iv) is to be replaced by the following condition at infinity: 
(v) The radiance remains bounded as the depth in the atmosphere tends to 
infinity 

lim J(r, fi, ip, no, ip ) < oo. (99) 

T — >00 

The main results for this model may be obtained from those derived for finite 
atmospheres by letting r b — > oo. All terms in Eqs. (12), (16) and (17) contain- 
ing F^(r b ) vanish, because they are multiplied by a factor that tends to as 
r b — > oo. Equation (21) thus becomes, in the semi-infinite case, 

S(r,fi ) = ^-B(T,fi ), (100) 

the function B for conservative semi-infinite media being given by (97). 
It is now possible to deduce the radiance at each point of the atmosphere. To 
that end, it suffices to substitute the present expression for the source func- 
tion in the formal solution (11)-(12) of the transfer equation. The analytical 
integration over r' then yields 



I(t, /i, ip, fio, y? ) = Y(-/i)F 5(ii + fJ, )5(ip - ipo) exp(r///) 

+ f H(, f^^ + ^^ , F . (101) 

This result has been known since the end of the fifties, as reported by Danielyan 
in a commentary to his Eqs. (5)-(6) in [28]. At the upper boundary r = 0, by 
(95) we retrieve the boundary condition (98) for — 1 < fi < 0, and also 

/(0.^o) = ^ g( ''> g( ''°V ,F. (102) 

47T fX + /i 

for < /i < 1. This is the well-known law of diffuse reflection for semi-infinite 
media [18], already retrieved by Stibbs [12] in the context of planetary atmo- 
spheres. It arises likewise from the finite-atmosphere counterpart (85) in the 
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limit of infinite optical thickness tj,, since the functions X and Y tend to H 
and 0, respectively, as T\> — > oo. As in the finite case, we might have established 
it directly, without resorting to the expression for the radiation field within 
the atmosphere. 

Deep within the atmosphere, the radiative field becomes, as indeed it should, 
asymptotically isotropic, with the limiting value obtained by letting r — > oo 
in (101) and by taking into account the fact that then 77 (r, fx) — > \/3. Conse- 
quently, 

Jim I(t, /J, if, /i , <p ) = —H(hq)hqF . (103) 

An isotropic radiance at infinity goes along with a vanishing net flux density at 
infinity, implying a vanishing net flux everywhere since the latter is constant. 
As a result, the upward and downward fluxes coincide everywhere, and their 
common value may be deduced from the expression (101) for the radiance. We 
get 



rl /id/i 

x / \ji v(t,Vo)+m(t,-v)}-^-- (104) 
Jo fj, + Ho 

For r = 0, the expression in brackets equals H(fi) by Eq. (95), and the so- 
called "alternative form" of the if -equation, as found for instance in [17], p. 
162, shows that the integral is equal to 2/H(/i ). One thus definitely retrieves 
F T (0,/i ) = F^O, /i ) = HoFq- In the deep layers of the atmosphere, as the 
function r\ tends to \/3, both fluxes approach the following limiting value: 

\/3 

lim F\t, /i ) = lim F l (r, /i ) = ^-H(fi )fi F . (105) 

The grayness factor of a semi-infinite atmosphere, defined as F^ (00, /i ) / liqFq 
by analogy with (61), is thus given by 

G(/i ) = ^H(no). (106) 

There exist numerous tables of the function H in the literature, for instance in 
the monograph by Chandrasekhar [18]. Our own tables reveal that G(/io) ex- 
ceeds unity for /j > 0.665, the maximum value 1.259 being reached for /x = 1. 
It may be noted that the expression (63) for the grayness factor of finite media 
does reduce to (106) as r b — > 00, since in that case a\ tends to the first-order 
moment of the function H, which has the value 2/\/3, while Pi approaches 
0. To our knowledge, expression (106) for the grayness factor was established 
for the first time by King [8] , who called it the "greenhouse factor" . We point 
out that in his derivation the factor /j is missing in both the numerator and 
denominator of the left-hand side of his Eq. (25), two omissions that happily 
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cancel each other. 

As to the temperature distribution of a semi-infinite atmosphere in local ther- 
modynamic equilibrium, it is readily found by simply identifying the source 
function S(t, /i ) in Eq. (100) with the spectrally integrated formula of Planck, 
i.e., (a/7r)T 4 (r, /i ). Remembering the relation crT 4 = F defining the temper- 
ature T , we arrive at the simple expression 



T(t, fi ) 



-B(t, ho) 



1/4 



(107) 



which is exactly the result that we would get if we let rj, — > oo in (79). The 
graphical representation of this temperature distribution is very much like 
that computed for a finite medium with r b = 10 (Fig. 1). As B(0, fx ) = H((j, ) 
and -B(oo,/i ) = VS/ioH (pLo) , the extreme temperature values are seen to be 
given by 



T(0,/xo) 
T 



1/4 



T 



^-/i H(ii ) 



1/4 



(108) 



from which follows the exact relationship T(oo, /i )/T(0, /i ) = (v^/A)) 1 ^ 4 - 
Equation (107) can be evaluated for grazing incidence of the incoming rays, re- 
minding ourselves that by virtue of the relations (95)- (97), the limit of B(r, /iq) 
for Hq — > + is H(0) = 1 if r = 0, and if r > 0. For grazing incidence, the 
temperature of the atmosphere is therefore zero everywhere except in the 
"layer" at r = 0, where it takes the value T /\^2. 

Now, for /i > 0, we may calculate the derivative of the temperature distribu- 
tion by substituting the expression (97) for B(t, /io) into (107), then appealing 
to Eq. (23) in [27], viz. 

Ho^(t, /i ) = 0(r) - 77(r, // ) (r > 0), (109) 

in order to get 

'^ = \H(^)[ ( f ) (r)-r ] (r^ )]. (110) 

The sign of T'(r, /i ) is therefore the same as that of the function <f>(r)— t)(t, /iq), 
that we have tabulated. We arrive at the following conclusion: 

- for < ii < 1, the temperature has a maximum in the interval ]0, +oo[, 
achieved at a depth r which is closer to the smaller /i , and then approaches 
T(oo,/i ) = [^/io^(/io)] 1/4 T , 

- for jiQ = 1, the temperature increases with r up to the value T(oo, 1) = 
[f#(l)] 1/4 T . 



T(t, Hq) 
To 
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8.3 Semi-infinite atmospheres with sources at infinity 



In the second category of semi-infinite models, suiting the giant planets, we 
choose a vanishing incident radiation as the upper boundary condition. We 
thus let 

(Hi) the incident radiation be naught at the top of the atmosphere 

J(0,/x) = (-l</i<0), (111) 
and choose "stellar-type" conditions at infinity: 

(v) in the deep layers of the atmosphere, the radiance I(t,li) increases more 
slowly than exp(r//i) for any /i > 0: 

hm I(r,n)exp(-r/n) = (0 < li < 1). (112) 

Thus we retrieve the model of a gray stellar atmosphere in radiative equi- 
librium that was studied by Hopf [24], Chandrasekhar [18] and others. The 
formal solution of the transfer equation reduces to the integrals in (11) and 

(12) , and hence the functions Jo, Fq and Fq on the right-hand side of Eqs. 

(13) -(15) vanish. We look for a solution, unbounded at infinity, of the homo- 
geneous integral equation (20) with u = oo. Hopf [24] showed that an infinite 
number of such solutions exist, all defined up to a multiplicative constant, but 
that only one corresponds to a given emergent flux F, viz. 

S(r) = —[r + q (r)]. (113) 

Here, q(r) denotes Hopf 's function, which may be expressed as [29] 

^^{ 1 + 5i'lM [1 - exp( - T/ ' ,)|d "}' (n4) 

One may then deduce the radiance at every point in the atmosphere, in much 
the same way as in Subsection 8.2. We get 



u \ 3F 
/(T ^ )= 4^ 



(115) 



The resulting radiance has been published often in the literature, but to our 
knowledge in forms more intricate than (115). At the surface of the atmo- 
sphere, Hopf's function has the value q(0) = l/VS, while 7/(0,//), given by 
(95), at once shows that 1(0, n) = for \x < 0. On the other hand, for fx > 0, 
we retrieve Hopf's result [24] 

I(0,n) = ^fHQi). (116) 
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As r — > oo, Hopf's function remains bounded, approaching q(oo) ~ 0.710, 
while the function 77 tends to \/3, so that 

3F 3F 
I(t,h) ~ — [r + g(oo) + //] ~ — r (r -> oo). (117) 

The radiance diverges at infinity, but without violating the condition (112). 
The upward and downward fluxes may be deduced from expression (115) for 
the radiance, much as in Subsection 8.2. We get 

3F \/3 r° 

F\r) = —[r + q(r)\ + \F J_^{t^)^ (118) 



F l (r) = ^[r + q(r)} - ^F V (r,^d^ (119) 

Note that the difference of these two fluxes just gives F, in agreement with 
(8), because, as may be shown, t](t, fi)fx 2 dfj, = 2/y/3. 
The calculation of the temperature distribution in an atmosphere in local ther- 
modynamic equilibrium, with the source function given by (113), is straight- 
forward; it leads to the famous relation of Hopf [24] 




(120) 



in which the effective temperature is defined by <rT t f > = F. 



9 Conclusion 



The generalized problem we have outlined originated with two prototype prob- 
lems of the classical theory of planetary atmospheres. The first one is of rel- 
evance to situations where shortwave radiation interacts with the planet's 
atmosphere, modeled as a slab under parallel-beam illumination from the sun 
and overlying a partly reflecting ground. In the second problem, there is no 
incident radiation on the upper surface of the atmosphere but only on the 
lower one, issuing from an isotropically emitting ground, a model appropriate 
to radiation transfer through a planetary atmosphere in the thermal (or in- 
frared) regime. In a strictly gray model, in which both the sun's radiation and 
the planet's ground emission are scattered and absorbed by the planet's atmo- 
sphere, the combined boundary conditions prevail simultaneously, providing 
us with the model that we offer in mathematical elaboration in the present 
article. 

Of course, this physically rather unsophisticated model is not realistic, but 
in exchange it is susceptible to analytical tractability. It yields the radiance 
field in a semi-infinite atmosphere illuminated from the outside or heated by 
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sources at infinity, for which the analytical developement could be pushed as 
far as it goes. In the case of an optically finite atmosphere, we were able to 
calculate exactly the source function for each depth, as well as the upward 
and downward fluxes across the boundaries. As a result, it was possible to 
determine the temperature of the atmosphere as well as that of the ground. 
Models more realistic than the one here presented may inspire further devel- 
opments and extensions. In particular, we think of the possibility of applying 
each set of boundary conditions to two separate frequency domains in such a 
manner that the radiation interactions remain coupled by the radiative equi- 
librium constraint. Models of that kind were proposed early, notably by Emden 
[3], and Milne [7], who realized that the angular variable no of the gray models 
could be extended to the range ]1, +00 [ in what is sometimes termed semi-gray 
models. This opens up the possibility of applying the results for gray media to 
those of semi-gray media, which represents a first step towards more realistic 
models. 

Semi-gray models are thus seen to constitute the natural next step in a hier- 
archy of analytical models of radiative transfer within planetary atmospheres. 
We intend to propose such a semi-gray model in two companion papers. The 
present study can therefore be judged as an intermediate step in the applica- 
tion of the exact solution to the thermal radiation problem, as given recently 
in [19], and the forthcoming semi-gray model. 



Appendix. Analytical expressions for the B- and £ _mnc tions in a 
conservative plane-parallel medium. 

The expressions below involve the auxiliary functions T, g and H defined by 
Eqs. (88), (90) and (86), respectively. We will also need the two basic functions 
(± = C±(aO already introduced for a finite slab in Section 6 of [30], as well 
as the moments of order and 1 of the X- and V-functions defined by (26). 
These moments are tabulated in [17], for instance. 
For u < 00 and \i > 0, we have 

B(t,h) = c{n)esxp(-T/n) + ^H(h)h[(_(-h)$ + (t, n) + ( + (-h)$-(t, 11)}, 

(Al) 

where c(/i) is given by (92) and 



+1 f 1 #r Uv)F*tlM. - exp( ~ (Tb " T)A V , (A2) 

2 Jo H(v) v-u v + a J K J 
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exp(-r/v) exp(-(r b - t)/v) 
v — jj, v + n 



}dv. 



(A3) 



These expressions simplify in a semi-infinite space, since then the functions 
3"± coincide with the ^-function defined by (93). Likewise, we have (± = 1 in 
a half-space, and Eq. (Al) reduces to Eq. (97). 

These formulae are new, they are given here without any proof. We obtained 
them by solving an integral equation of Cauchy type, satisfied by the function 
B, along the general lines of [31]. Our equation (Al) is a particular case of 
the general equations in Section IV of [31], from which it follows by choosing 
a = 1 and cq{z) = exp(— r/z). 

The expression of the £ -function has already been given in [19], Eqs. (62)- (63), 
which we repeat here for convenience. For rt < oo 



while for 77, = oo, we have £o( r ) = 1 as expected. 
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